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ABSTRACT 



u 

The Horava - Lifshitz (HL) theory has recently attracted a lot of interest as a viable 
Q . solution to some quantum gravity related problems and the presence of an effective 

' cosmological constant able to drive the cosmic speed up. We show here that, in the 

C/2 | weak field limit, the HL proposal leads to a modification of the gravitational potential 

because of two additive terms (scaling respectively as r 2 and r~ 4 ) to the Newtonian 
1/r potential. We then derive a general expression to compute the rotation curve of 
an extended system under the assumption that the mass density only depends on the 
■ cylindrical coordinates (R, z) showing that the HL modification induces a dependence 

■^j- | of the circular velocity on the mass function which is a new feature of the theory. As 

a first exploratory analysis, we then try fitting the Milky Way rotation curve using its 
visible components only in order to see whether the HL modified potential can be an 
alternative to the dark matter framework. This turns out not to be the case so that 
(y-) . we argue that dark matter is still needed, but the amount of dark matter and the dark 

C^> halo density profile have to be revised according to the new HL potential. 



o 
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1 INTRODUCTION der to reduce the number of independent terms entering the 

action, the principle of detailed balance was used. Unfortu- 
Inspired by the Lifshitz theory in condensed matter physics, a goon became ^ that tMg second condition leads 

Horava (2009a,b) has recently proposed a new theory of tQ problems in the low energy limit (Lu et al . 200 9) thus 

gravity based on an anisotropic scaling of space and time .. ,. , c A . c c • ■ i UT 

b J o±* motivating the search lor modification of the original fiL 

in the UV limit. Usually referred to as Horava - Lifshitz theQry where the actiQn breakg ^ detailed b£jance CQndi _ 

(hereafter, HL) theory, the HL proposal shows a reduced tiQn either goftly (Kehagias & sfestos 2009; Lee e t al. 2009; 

invariance, dubbed foliation preservzng dzffeomorphism in- Capasgo & Polychronakos 2009; Kiristis & Kofinas 2009a; 

variance, which however reduces to the standard one in the Kiristig & Kofinas 200gb) Qr nQt In particular> SotirioUi 

IR limit where General Relativity is recovered. As an attrac- yisser and Weinmrtner (2009a)b) have worked out a modi _ 

tive feature, the HL theory turns out to be power counting c , UT .,, j , ., j u , 

' J fa nec j theory with no detaiied bafance condition imposing 

renormalizable which has motivated the great interest in in- , „ n • , , , i 

& that oniy parity preserving operators enter the potentiaf. 

vestigating with great detail its theoretical and cosmological 

aspects. It is worth remembering that, in its original formu- Although the discussion about the foundations and the 
lation, two conditions were imposed in order to drive the possible conceptual and phenomenological problems of the 
choice of the field action. First, the projectability condition HL theory and its modified versions is still open, it is never- 
was supposed to hold true. Since time plays a fundamental theless worth systematically investigating its consequences 
role from the very beginning, it was assumed that the lapse at every scale. In particular, it is interesting to study its 
function (defined when working in the ADM formulation of static spherically symmetric solutions since one can thus de- 
gravity) has to be a projectable function on the spacetime rive the gravitational potential generated from a point mass 
foliation, that is to say a function of time only. Second, in or- source. Recently, this problem has been addressed by Tang 
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& Chen (2009) for the HL theory with the projectability 
condition and no detailed balance. They argue that only 
the Minkowski or de Sitter spacetime are solutions, but wc 
will show here that this is actually not the case. As a conse- 
quence, we find that the gravitational potential generated by 
a point mass differs from the Newtonian one because of the 
presence of additional terms depending on the HL coupling 
parameters. Such terms have to be taken into account when 
computing the potential generated by an extended system, 
such as a galaxy. Therefore, we work out a general formal- 
ism to estimate the rotation curve (i.e., the circular velocity 
v c as function of the distance R from the centre) for an ex- 
tended source showing that the HL theory can boost the 
v c (R) with respect to the Newtonian value. Motivated by 
this consideration, we then try to fit the Milky Way rota- 
tion curve using visible matter only to both constrain the 
HL parameters and investigating whether it can work as an 
effective dark matter component. 

The plan of the paper is as follows. In Sect. 2, we look 
for static spherically symmetric solutions for the HL theory 
with projecatibility condition and show that there is indeed 
a new solution leading to a modified gravitational poten- 
tial. A general formalism to compute the rotation curve for 
an extended system is presented in Sect. 3 and then used 
in Sect. 4 to work out the predicted rotation curve for the 
Milky Way. Here, we also present the data and the results of 
fitting them with our modified potential and no dark matter. 
Conclusions are finally given in Sect. 5. 



2 THE POINT MASS POTENTIAL 

The weak field limit of the HL theory as proposed by 
Sotiriou, Visser and Weinfurtner has been yet discussed by 
Tang & Chen (2009, hereafter TC09) so that here we will 
only summarize the main steps and stress where our work 
differs from their one. As stressed in TC09, we start from 
the observation that not all the static spherically symmetric 
solutions found for (modifications of) HL theory preserve 
the projectability condition. In fact, if one assumes a metric 



ds 2 = -N 2 (r)dt 2 + +rW 
9(r) 



(1) 



1, or simply reduces to Minkowski by means of a redefinition 
of the radial coordinate. The second choice seems therefore 
more interesting. In such a case, the equations obtained by 
varying with respect to / and N reduce to 



dN 2 N 2 N 2 (t) , . 



dr 



Jo 
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+ - r/' (4g4 + 2ff 5 +ff6) 



+ 



8/(1 -ff 



(97 + 98 ) 
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For / = 1, we get vi(r) = v-zir) = —go(, 6 r so that Eqs.(4) 
and (5) are equal hence the solution for N r (r) is the same. 
However, differently from what stated in TC09, this is not 
the only possibility. Indeed, in order to have the same N r (r) 
solving both Eqs.(4) and (5), one must have vi(r) = i^(r) 
which is possible by equating the coefficients of the terms 
with equal orders in r in the two functions. Comparing v\ (r) 
and v%(r), one has just to equate the coefficients of the terms 
in r~ 3 and r -5 thus obtaining the following two equations : 

(l-/)[2 S2 (l + 7/)+ S3 (l + 5/)] = (W) 2 (2 52 +4 53 ) ,(8) 



mation dt = dtpa — [(1 — N 2 ) 1 '" 2 /N 2 ]dr, one can show that a 
necessary condition for the projectability condition to hold 
is that g = N 2 . The line element therefore reads : 



-N\t)dt 2 + 



2 (dr + N r dt) (dr + N r dt) 



f 



+ r'dQ. 2 



(2) 



The equations of motion for such a metric are given in 
TC09 and will not be repeated here for sake of shortness. 
The authors then argue that in the IR limit the / func- 
tion is constrained to be 1 thus ending up with the usual 
Schwarzschild - de Sitter metric as only solution. We here 
show that this is actually not the case. To this aim, we first 
note that the Euler - Lagrange equation for N r — N r /f in 
the IR limit reduces to 

f_Nr 

f r 

which is solved by either N r — or / = const. The first 
choice gives Minkowski back when choosing go = and / = 







(3) 



2(1 - /) 2 [4 ff4 (l + 23/) + 2 55 (1 + 17/) + pb(1 + 14/)] + 
8/(1-/) [2^(1 + 7/) +9„(1 -4/)] = 
2(1 - /) 3 (4 3 4 + 2g 5 + go) + 8/(1 - ff (g 7 + g 8 ) , (9) 

which allows us to set two of the quantities (gi, . . . , gs, /) as 
a function of the others. Provided this condition has been 
satisfied, the solution of Eq.(4) automatically solves also the 
integral condition (5) so that we can limit our attention only 
to Eq.(4). This is a linear first order equation in N 2 (r) which 
can be analitically solved giving : 



*H-W3-£(f« + B-°-J> 



(10) 



r 2 2p V3 r 
where M is an integration constant and we have defined* 



* Note that, following TC09, we have used units in which Z = 1 
with Z a dimensional parameter of the HL theory. In order to 
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( A = ~ 9 oC b 

B = 2C 4 (1-/) 

C = 2C 2 (l-/)[2<7 2 (l + 7/)+ 5 3(l + 5/)] 

D = 2(l-/) 2 [4 54 (l + 23/) + 2p 5 (l + 17/) 

+ ff6 (l + 14/)]+8(l-/)/[2< ?7 (l + 7/) 



(11) 



+ 38(1-4/)] 

In order to derive the gravitational potential, we have first 
to rewrite the line element in the usual Schwartzschild - like 
form, i.e. 

ds 2 = — goo{r)dt 2 + g rr (r)dr 2 + r 2 dQ 2 



and then use the general relations : 
ff oo(r) = N 2 - fN 2 , g rr (r) 



N 2 



f(W-fN?) ■ 

Fixing N = 1 (which is always possible by a rescaling of the 
t coordinate), we then get : 

'A 

r If {V 
g rr (r) = l//flooM . 



9oo(r) = 1 - tl± + 1 ( + U 



_C_ _ D_\ 

3r 4 J ' 



(12) 
(13) 



Since the gravitational potential generated by a pointlike 
mass particle may be easily recovered from the usual relation 
5oo(t*) = 1 + 2$(r), we then find : 



(14) 



Note that, up to now, the two integration constants / and 
M are still undetermined. However, in order to recover the 
usual Newtonian potential in the GR limit (i.e., for A = B = 
C = D — 0), we must set fM = 2Gm with m the mass of 
the gravitational field source. It is worth noting that, for 
m = 0, the potential (14) leads to unphysical divergences 
in the metric because of the terms in C/r~ 2 and D/3r 4 . 
In order to avoid this problem, we must impose that both 
C and D vanish for m = 0. The first trivial choice is to 
impose C — D — identically so that one recovers the 
usual Schwartzschild - de Sitter case. As a more attractive 
possibility, one can postulate a dependence of the metric 
coefficient / on the mass. In such a case, Eqs.(8) and (9) 
can not be read as a relation among / and the HL couplings 
<7i since this will induce a dependence of the HL Lagrangian 
on the mass which is not possible. As such, we consider them 
as algebraic equations for / so that, in order to be fulfilled, 
we must equate the terms with equal powers of / on the two 
sides. We then obtain : 

92 = gz = 

<? 4 = -(36<? 5 + 15 36 + 4(? 7 )/96 
t ff8 = (16/3) fl7 



have (A, B, C, D) expressed in the more common c = 1 units, 
one has simply to multiply each term in Eq.(ll) by with Q 
having, in these units, the dimension of a length. 



which implies C — and 
(I-/) 2 



D 



12 



[12(1 - f)g 5 - 9(1 - f)g 6 + 4(1 - 153/) fl7 ] 



With these conditions, we can finally write the gravitational 
potential generated by a points mass m as : 



$(r) = <5jv(r) + <S>HL{r) 



(15) 



where <I>jv(r) = —Gm/r is the Newtonian potential and 

A 2 B D 1 
*««=12/ r + -f-m~ 4 (16) 

is the correction due to HL theory. Is is worth noting that, 
while the second term is simply a constant having no impact 
on the dynamics, the other two terms have a simple relation 
with previous results in literature. Indeed, the first one plays 
the role of a cosmological constant thus remembering the 
Schwartzschild - de Sitter solution already found in TC09. 
On the other hand, the third term has the same asymptoti- 
cal behaviour of the corrections to the Newtonian potential 
in the Kehagias - Sfestos (2009, hereafter KS) static spher- 
ically - symmetric solution of HL gravity. Iorio & Ruggiero 
(2009) have shown that such corrections are proportional 
to m r~ thus suggesting that D should actually be a func- 
tion of the mass m generating the gravitational field, a point 
which we will come back to later. 

The corrective term may be conveniently written as : 



GMr. 



Br s 



4/GM 



+ 



1JD 



(17) 



with M© the Sun mass and r s an arbitrary chosen reference 
radius introduced to define the dimensionless quantity r) — 
r/r s . In Eq.(17), we have finally defined the scaling radii : 



I At s ) 



Ar s 

Dr~ 3 
12/GM G 



1/4 



(18) 



which can be used as model parameters instead of the (A, D) 
coefficients. Some caveats are in order here, First, note 
that, in order to preserve the interpretation of (rA,ro) = 
r s x (r}A,r)D) as physical radii, one has to postulate that 
the (A, D) parameters are positive quantities thus narrow- 
ing the space of the parameters (go, 55, <?6, gi,f) entering 
(A, D) . Should A or D be negative, we could nonetheless de- 
fine (rA,rc) as above and accordingly change the sign of the 
corresponding term in the potential. For definiteness, we will 
both A and D are positive so that the potential is given by 
Eq.(17) without any sign change. We then remember that, 
in our scheme, / is a dimensionless function of fi — m/M®, 
but its functional dependence can not be obtained in any 
way. The only constraint we have is f(fi = 0) = 1 so that 
D = and we recover the Schwartzschild - de Sitter solution 
in accordance with the fact that go (and hence A) play the 
role of a cosmological constant term. In order to parametrize 
our ignorance we can redefine the above scaling radii as : 
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(\2f&GMo\ 1/2 
\ Ar s ) 



1 1/2 



/© 



T)D = 



Dp,r: 



12/qGMq 



i" 




/© ' 


1/4 




_ D Q 


/(/*)_ 





(19) 



where quantities labelled with a are evaluated for p — 1. 
It is worth wondering whether some hint on the functional 
expression of f(p.) can be retrieved. Up to now, we have 
postulated / = f(p) in order to get a mathematically viable 
solution other than the Schwartzschild - de Sitter one. It is, 
however, worth noticing that a dependence of / on the mass 
may also be physically motivated. Although the HL theory 
is obtained by modifying General Relativity, it is still true 
that the properties of the spacetime are determined by the 
source of the gravitational field. As such, one can expect that 
the corrective HL term to the Newtonian potential is still 
related to the only property characterizing the field, i.e. the 
source mass m. Since the (A,B,D) coefficients in Eq.(??) 
are related to the HL couplings (and hence are universal 
quantities), the only way to introduce a dependence of the 
solution on the source properties is to postulate / = f(p). 
Although qualitative, this discussion shows that our math- 
ematical assumption is actually deeply related to a physical 
motivation. 

Finally, we note that the second term in Eq.(17) simply 
adds a constant to the potential which has no effect in any 
situation of interest so that we will henceforth neglect this 
term. Note that this by no way means that B can be set 
to zero. Indeed, B = means / = 1 so that also D vanish 
and we go back to the Schwartzschild - de Sitter solution, 
while we are here interested in the more general case. We 
therefore assume B 0, but nevertheless neglect its contri- 
bution hereafter because drops off from the derivation of the 
quantities we are interested in. 



3 THE ROTATION CURVE 

The modified gravitational potential derived above devi- 
ates from the Newtonian one because of the additive terms 
in Eq.(17). Depending on the values of the scaling radii 
(rA,rD), we can have different situations. As a general re- 
mark, we note that, for r\ « 1, the last term in Eq.(17) 
increases the potential with respect to the Newtonian one. 
On the contrary, for r\ » 1, the first term boosts the poten- 
tial making it deviating from the Keplerian fall off. When 
considering the rotation curve, v c (r) = rd<&/dr, we thus get 
a circular velocity which may significantly differ from the 
Newtonian one being larger than the classical value both in 
the inner and outer regions. It is therefore worth wondering 
whether such deviations may help in fitting spiral galaxies 
rotation curves without the need of other mass than the 
visible one. That is to say, we are here interested in inves- 
tigating whether the HL modifications to the potential can 
also play the role of an effective dark halo. As a preliminary 
remark, we stress that a similar analysis could in principle 
be made also for the KS solution. However, in that case, the 
corrections fade away as r" 4 and there is no cosmological 
constant term. It is therefore expected that the correction 
to the rotation curve in the outer regions are negligible so 



that we argue that the KS solution can not play the role of 
an effective dark halo. 

To this end, we have first to derive an expression for the 
rotation curve of an extended system generalizing the proce- 
dure adopted in the Newtonian gravity framework. In that 
case, the circular velocity in the equatorial plane is given by 
Vc(R) = Rd$>/dR\ z= o, with $ the total gravitational poten- 
tial. Thanks to the superposition principle and the linearity 
of the point mass potential on the mass m, this latter is com- 
puted by adding the contribution from infinitesimally small 
mass elements and then transforming the sum into an inte- 
gral over the mass distribution. Such a simple procedure can 
not be applied in the HL case since the corrective term $ai 
does depend on m in a way that we do not explicitly know 
so that a nonlinear dependence can not be excluded a priori. 
In order to overcome this difficulty, we have developed an 
alternative procedure which can be actually applied to any 
kind of potential provided some general conditions hold. 

As a starting point, let us denote by F p (m, r) the grav- 
itational force (per unit of test mass particle) generated by 
a point mass m. Whatever is the dependence of F p on m, it 
is always true that the total force due to N particles of mass 
m is the sum of the single forces. As such, taking the con- 
tinuum limit, we can estimate the (magnitude of the) total 
force as : 



F(r) = J n(m, r')F p (m, |r - r' 



\)dmdV 



(20) 



where the integral is over the full mass range and volume 
and n(m,r) is the star mass function (hereafter MF), i.e. 
the number of stars in the volume element dV with mass 
between m and m + dm. Because of its definition, we have : 



J m. 



n(m, r)mdm — p(r) 



(21) 



with p(r) the mass density. As usual in literature, we will 
adopt the factorization hypothesis thus writing n(m, r) = 

ip(m)p(r) with ip(m) the local ^ MF and p = p(r)/po with 
po = p(Ro). Defining p = m/M , we have the following 
normalization condition for the local MF : 

pip{p)dp = po/M Q 

so that henceforth we use as local MF the quantity J\fiip(p) 
with 

l 

(22) 



J y-min 



Ml 

Let us now assume that the point mass gravitational force 
may be factorized as : 

F P (p,r) = ^Mp)Mv) 

>s 

with fn and f r {r) dimensionless functions depending on the 
particular form of the point mass gravitational potential $ p . 
Remembering that F p = — V<I>j,, it is only a matter of alge- 
bra to show that it is : 



T Note that the term local refers to the Solar neighborhood when 
the galaxy is the Milky Way. For external galaxies, by local, we 
mean the MF in the neighborhood of a suitably chosen reference 
radius Rg. 
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= m > fr(v) = i/v 2 , 

for the Newtonian potential. For the HL corrective potential, 
we can split it as the sum of two terms, i.e. 

®HL(r,r A ,r D ) = $A(r,r A ) + ^^(r,^) , 

and then obtain : 

7(m) 



'•-(-)(-*-") 



for $a and 

/r = 



/0 

^0 /(m) 



for the $£> term. We now use cylindrical coordinates (7?, 9, z) 
and the corresponding dimensionless variables (rj, 9, (with 
£ = z/r s ) and rely on the factorization hypotheses for both 
the MF and the point mass force to finally get : 

F(r) = G po r s 

poo POO />7T 

x / v 'dr,' / dc' / MA)p( v , ,e',c')do' (23) 

J J —00 Jo 



21 1/2 



(24) 



with the shorthand notation 

A = + v ' 2 - 2W cos (61 - 9') + (C - O 

The circular velocity in the equatorial plane along the ma- 
jor axis (which is the quantity typically measured for spiral 
galaxies ) will be simply v 2 (R) = RF(R,9 = z = 0). Since 
we will be interested in axisymmetric systems, we can set 
p — p{rj,Q. Moreover, our systems will be spiral galaxies, 
hence made out of a spheroidal bulge and a circular disk, so 
that a convenient choice for the scaling radius r s will be the 
disk scalelenght Rd- Under these assumptions, the rotation 
curve may then be evaluated as : 



v 2 c (R) 



G po R 2 , 7} 

/ V'dv' / P(v',(')d(' 

J0 J - oo 



/r(A ) 



(25) 



with 



121 1/2 



(26) 



A = A(9 = C = 0) = [r, 2 + T)' 2 - 2r)r)' cos9' + £ 

It is worth noting that, except in very particular cases (e.g., 
the Newtonian potential for a spherically symmetric mass 
distribution), the integrals over the coordinates in Eq.(25) 
have to be evaluated numerically. For computational rea- 
sons, it is useful to exchange the order of integration and 
resort to the logarithmic variables A = log 77 and cj = log£ 
so that we get the following equivalent expression for the 
rotation curve : 



v 2 (R) = G po R 2 d (lnlO) 2 dex(A) 



with dex(;r) = 1Q X and 



d(, 



(27) 



dex(2A')dA' / f r (A )p(X', w')dex(w')dw' . (28) 



R = 



Eqs.(25) and (27) are fully general and can be used to com- 
pute the rotation curve provided the expression for /^(/j) 
and f r (ri) are given. As a consistency check, it is easy to 
show that, for the Newtonian potential, the term depending 
on the MF is identically unity so that Eq.(25) reduces to a 
simple rewriting of the standard result. For the HL term, we 
get a dependence on the MF through the multiplicative term 
on the second row of Eq.(25). It is worth stressing that the 
MF here only plays the role of scaling up or down the rota- 
tion curve. This is a consequence of the HL corrective term 
(17) not depending on m. As such, indeed, one has simply to 
sum the contribution of all the stars notwithstanding their 
mass and this is indeed what the MF term gives in the HL 
case. The total rotation curve for the HL theory will finally 
be given as : 

v 2 (R) = v 2 N (R) + v 2 HL (R) = v 2 N (R) + v\{R) + v 2 D {R) 

where one can use the literature results for the computation 
of the Newtonian term v% (R) and our general rule with / M = 
1 and f r (v) given above to estimate the HL contribution 
vjj L (R). Note that, in order to evaluate v 2 (R), we need to 
know not only the mass density (as in the Newtonian only 
case), but also the local MF because of the HL term. 



4 THE MILKY WAY ROTATION CURVE 

The HL theory has been originally conceived as an attempt 
to solve some quantum gravity related problems, but has 
soon attracted a lot of interest even as an alternative model 
for dark energy because of the presence of a cosmological 
constant like term. We have here shown that also the grav- 
itational potential is modified with respect to the standard 
Newtonian one so that it is worth wondering if the addi- 
tional terms may help in reconciling the data on the rotation 
curves with what is predicted from the visible matter only. 
To this end, one should fit the rotation curves of many spi- 
ral galaxies and find out that the fit is indeed satisfactorily 
good and the model parameters (ta, td) are the same for all 
the galaxies being related to the HL Lagrangian couplings 
and hence universal quantities. Such a task is actually quite 
complicated because of the need for the knowledge of the lo- 
cal MF which is not at all constrained in external galaxies. 
Moreover, what one observes for these systems is the sur- 
face brightness profile due to the two visible components, 
namely the bulge and the disk. What we need to evaluate 
the rotation curve is, however, the mass density so that, 
even assuming simple functional forms for the bulge and 
the disk, we still do not have any knowledge of the bulge 
and disk mass - to - light ratios thus adding more parame- 
ters (and hence severe degeneracies) to be determined. As 
a first preliminary test, we therefore limit our attention to 
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the Milky Way (hereafter MW) only. Our position within 
it allows us to determine the local MF (e.g., by star counts 
or converting the observed luminosity function into a MF 
through an empirically determined M/L ratio) thus reduc- 
ing the uncertainties of the problem. Moreover, we also have 
direct determinations of the galactic parameters of interest 
so that the only unknown quantities are the HL parameters 
{rA,ro) thus strongly reducing the possibility of degenera- 
cies. In the following, we first describe the MW mass models 
and the data on the rotation curve and then present the fit- 
ting procedure and the results. 



4.1 The mass models 

It is common to describe a spiral galaxies as the sum of a 
visible component (made out of stars and gas) embedded in 
a dark halo (mainly populated by cold dark matter parti- 
cles). Lacking any definitive laboratory evidence for DM, the 
only reason why one has to include it in galaxy modelling 
is to fit the rotation curves data. Since this evidence im- 
plicitly assumes that the Newtonian potential theory is the 
correct one, there is actually no compelling reason why one 
should add a priori a dark halo to the visible components in 
a modified gravity framework as the current HL theory. We 
therefore model the MW as made out of vibile matter only 
distributed in a spheroidal bulged and a thick disc. 

We follow Dehnen & Binney (1998, hereafter DB98) 
describing the bulge as a truncated power -law model, i.e. 
the scaled mass density reads : 



*<*-> = (£H 1+ l; 



7-/3 



CX P ~T72 



(29) 



with g 2 = R 2 + z 2 /q 2 . Fitting the model to the infared pho- 
tometric COBE I DIRBE data yields values for the model 
parameters, namely : 



= 7 = 1. 



0.6 



Rb = 1 kpc , R t = 1.9 kpc 



The reference density is not determined from the photome- 
try, but may be related to the total bulge mass Mb as : 



Po,b 



Mb 



2nf °°RdRj~ oo p(R,z)dz 



(30) 



Dwek et al. (1995) has found M b = (1.3 ± 0.3) x 10 10 M 
for the bulge mass so that we set M b to its central value 
neglecting the measurement uncertainty. We have then to 
set the local MF for the bulge. Zoccali et al. (2000) have 
determined it from the luminosity function of lower main 
sequence stars finding : 

V>&(/i) oc p _/3 

with /3 = —1.33 ± 0.07 for stars in the mass range 
(0.15,1.0) M . We extend it to the brown dwarfs re- 
gion and neglect the measurement uncertainty thus setting 



t Actually, there arc different evidences that the bulge has triax- 
ial structure. We nevertheless use a less detailed spheroidal model 
since the bulge contribute to the dynamics over the range probed 
by the data is much smaller than the disc one, independently on 
the gravitational theory adopted. Such a simplification allows us 
to use Eq.(25) without introducing any significative bias. 



(P,fimin,fJ.max) = (-1.33,0.03,1.0) as our bulge MF param- 
eters. 

While important in the inner regions, the bulge only 
plays a minor role in determining the circular velocity over 
the regions probed by the data. The dynamics is here dom- 
inated by the disc component that we model as a double 
exponential. The mass density then reads : 

M*,*)=exp(-|-) e xp(-M) (31) 

where we follow DB98 fixing : 
Rd = KdRo , z d = 0.18 kpc , 

with K d — 0.30 ± 0.05 a scaling parameter and Ro the Sun 
distance from the MW centre. We then follow Cardone & 
Sereno (2005) adopting n d = 0.30 and Ro = 8.5 kpc as 
fiducial parameters. The disk reference density is related to 
the disk total mass as : 



Po,d = 



M d 



^R 2 z d 

while the disk mass is estimated as 



(32) 



M d = 2tt^E exp (R /R d ) (33) 

with E Q = 48 ± 8 M Q /pc 2 (Kuijken & Gilmore 1989) the 
disc surface density at the Sun position. A caveat is in order 
here. The measured value of E© refers to the total mass, 
both stars and gas. In principle, we should separate the two 
components computing the total disc mass using E* instead 
of E Q , with E* = E Q - E/sm and Ej SM = 14 -5 M /pc 2 
(Oiling & Merrifield 2001) the gas surface density. We should 
then add a further disc like component for the gas with 
its density profile. As a reasonable approximation, we can, 
however, assume that the gas follows the same density profile 
as the stellar disc with the same scalelength radius so that 
the total disc mass is indeed given by Eq.(33). 

As a final ingredient, we need the local disc MF. This 
is determined quite accurately thanks to our position in the 
disc plane. We refer the reader to Chabrier (2003) and refer- 
ences therein for a detailed discussion of this issue motivat- 
ing our choice of adopting the Kroupa MF, i.e. a power - law 
MF with slope j3 changing with the mass range as : 

0.3 0.01 < n < 0.08 



= < 



1.3 0.08 < fi < 0.50 
2.3 0.50 < n < 1.0 

I 4.5 1.0 < /x < 120.0 

The four different range contributes to the local mass density 
according to the following percentages : 

{h,h,h,U) = (7.39,48.21,29.22,15.18)% 

so that we can compute the multiplicative term entering the 
rotation curve for the HL term. 



4.2 The data 

Although not ideal targets for our test because of the poor 
knowledge of their parameter and MF, external galaxies are 
better suited for the measurement of precise and extended 
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rotation curves. On the contrary, our position in the MW 
disc equatorial plane makes it difficult to observationally 
determine this quantity. Indeed, one has to characterize the 
full line of sight velocity distribution in order to correct the 
observed velocity for asymmetric drift and projection effects. 
Moreover, the distance to the tracer should be known with 
great accuracy not to bias in a dangerous way the estimated 
circular velocity. Notwithstanding these difficulties, the ro- 
tation curve in the outer regions have been measured relying 
on Cepheids and HII regions. We follow DB98 to estimate 
the circular velocity from the Cepheids data of Pont et al. 
(1997) and the HII molecular clouds sample of Brand & 
Blitz (1993) using their same selection criteria. This dataset 
probes the radial range 8.2 < R (kpc) < 18.94 and is af- 
fected by a large scatter due to both measurement errors and 
the inhomogeneity of the data. In order to smooth the data 
without introducing any bias or spurious correlations, we use 
the local regression method (Loader 1999). Originally pro- 
posed by Cleveland (1979) and further developed by Cleve- 
land and Devlin (1988), the local regression technique com- 
bines much of the simplicity of linear least squares regression 
with the flexibility of nonlinear regression. The basic idea re- 
lies on fitting simple models to localized subsets of the data 
to build up a function that describes the deterministic part 
of the variation in the data, point by point. Actually, one is 
not required to specify a global function of any form to fit a 
model to the data so that there is no ambiguity in the choice 
of the interpolating function. Indeed, at each point, a low de- 
gree polynomial is fit to a subset of the data containing only 
those points which are nearest to the point whose response is 
being estimated. The polynomial is fit using weighted least 
squares with a weight function which quickly decreases with 
the distance from the point where the model has to be re- 
covered. We hence use this method § to smooth the sample 
points and cut data with R > 14 kpc since the local regres- 
sion method becomes unreliable for these points because of 
the sparseness of the sample in this region. We finally give 
it away points with S/N < (S/N) m in where the threshold 
signal -to -noise ratio has been set to {S/N) m i n ~ 8 for rea- 
sons explained below. We will refer to this sample as the DB 
data in the following. 

Another possible tracer of the velocity field in the outer 
regions is represented by the Blue Horizontal Branch (BHB) 
stars. Photometry and spectra from the Sloan Digital Sky 
Survey (SDSS) for a sample of 2466 BHB stars allow to de- 
termine the rotation curve up to ~ 60 kpc as described in 
Xue et al. (2008, hereafter X08). In order to extract v c (R) 
taking care of all the possible projection effects and sys- 
tematic biases (due to, e.g., the survey mask and targetting 
algorithm), X08 have relied on mock matched observations 
based on two different cosmological simulations of MW-like 
galaxies. As such, there are two different datasets labelled 
as N and S depending on the simulation adopted. We here 
use only the N dataset since it is in better agreement with 
the DB sample over the range of overlap. Note that this 
sample (referred to hereafter as the SDSS dataset) is made 
out of only 10 points obtained by binning the 2466 BHB 
stars in almost equally spaced radial bins covering the range 
(7.5,55) kpc. The S/N is quite high with a median value 



See, e.g., Capozzicllo et al. 2007 for a step - by - step description. 



{S/N) med — 8 so that we set (S/N) min = (S/N) med in or- 
der to give the same statistical weight to both datasets. As 
a final remark, we note that, while the DB data probe only 
the outer disc being R/Rd ^ 3.2 — 5.4, the SDSS sample 
extends mainly in the halo region most of the data having 
R/Rd > 4 and extending up to R/Rd — 20. As such, the two 
sample nicely complement each other and allow us to check 
whether the outer rotation curve may be reproduced with- 
out any dark matter contribution and, at the same time, still 
preserving the agreement with the data in the disc region. 



4.3 Fitting procedure 

In order to constrain the HL model parameters, we employ 
a standard Bayesian approach first defining the likelihood 
function as : 



£(p) tx exp 



= exp ■ 



x 2 (p) 



f s (Ri) 



(Ri 



(34) 



where p = (logr/A, log 77c, log 770) are he HL model 
parameters^, v° bs (Ri) and v t J l {R i ) are the observed (with a 
measurement error e/) and theoretically predicted values of 
the circular velocity at the radius Ri of the i - th point and 
the sum runs over the 101 DB and 10 SDSS data points. 

The best fit is obtained by maximizing the likeli- 
hood £(p), but it is worth stressing that, according to the 
Bayesian philosophy, the best estimate of the parameter pi 
is not the best fit one. On the contrary, one has to marginal- 
ize over the remaining parameters and look at the shape of 
the marginalized likelihood function defined as : 

Ci(pi) oc J £(p)dpi . . . dpi-idpi+i . . . dp n 

with n the total number of parameters. Actually, what 
we do is running a Monte Carlo Markov Chain code to 
efficiently explore the three dimensional parameter space 
(log rjA, log 77c, log rjo) and use of the histogram of the val- 
ues for the parameter pi to estimate the mean, the median 
and the 68 and 95% confidence ranges. Note that, because of 
degeneracies among the model parameters, the best fit pa- 
rameters pi,/ may also differ from the maximum likelihood 
ones pml, i.e. the set obtained by maximizing each of the 
marginalized likelihood functions. 



4.4 Results 

As a preliminary discussion, it is worth clearly stating which 
are the parameters we can constrain. Eq.(25 shows that, 
in order to compute the HL contributions to the rotation 



^ We use logarithmic units in order to explore a wider range. 
Note also that, rigorously speaking, Eq.(34) is not correct since 
the SDSS data points are somewhat correlated being obtained by 
a binning procedure. However, since the bin spacing is quite large, 
it is likely that the correlation matrix (not available) is close to 
diagonal so that Eq.(34) is essentially correct. 
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Figure 1. Best fit HL model rotation curve superimposed to the 
DB and SDSS data. 

curve, namely v A {R) and Vp(R), we must know the func- 
tional expression of f(fi) entering /^(/i) through the scaling 
radii (va^d)- Since we do not know this function, we can 
not separately constrain the parameters (r^©, tdq)- It is, 
however, easy to show that : 

kav a (R,ija,@) = v A (R,r/Ai) , 

K D vl,(R,r)D,e) = vl,(R,r]Di) , 

where we have defined : 

while rjA = rA/Reff and r/u = ro/Reff are considered as 
constant parameter not depending on fi and we assume that 
v A (R,riA) and v 2 D (R,r)D) are evaluated setting / M = /i in 
Eq.(25). Our fit will then give constraints on (tja^d), while 
those on (VAEnVDo) could be derived provided a theoreti- 
cally motivated functional expression for f(j-i) is given. 

Skipping to logarithmic units to investigate a larger 
range, we run a single chain with 150000 points reduced 
to ~ 10000 after cutting out the initial burn in phase and 
thinning the chain (with a step of 10) to avoid spurious cor- 
relations. The best fit point turns out to be : 

(log VA , log r) D ) = (1.051,-0.060) 

giving x 2 /d-o.f- = 523.56/109 = 4.80. Needless to say, such 
a high reduced x 2 is a strong evidence against the HL model 
with no dark matter. A simple look at Fig. 1 shows that 
the large \ 2 is due mainly to the SDSS data. Indeed, it 
is X 2 D B /d.o..f. = 235.91/99 = 2.38 and XsDSs/d.o.f. = 
287.65/8 = 35.96 for the DB and SDSS samples, respec- 
tively. Although the model is clearly at odds with the data, 
we give below, for completeness, the constraints on the pa- 
rameters from the analysis of their histograms along the 
chain. For log tja we get : 



(log T) A ) = 1-049 , (log ?M)med = 1.050 , 

68% : (1.040, 1.059) , 95% : (1.028, 1.066) , , 
while it is : 

(log VD ) = -0.061 , (log j) D ) m .d = -0.060 , 

68% : (-0.064, -0.056) , 95% : (-0.068, -0.054) , 

for logrjo- Converting median values in linear units, we get : 

{r A ,r D ) = (11.2,0.87) kpc 

which allows us to make some qualitative interpretation of 
the results. Considering the form of the potential for the 
point mass case, we see that only the first term can actually 
boost the stellar contribution to the rotation curve. This 
term essentially translates in a contribution to v c (R) which 
linearly increases with R/va so that ta must be small in 
order to fit the data in the DB probed region without dark 
matter. On the other hand, a too small value may lead to a 
v c (R) soon diverging thus spoiling down the agreement with 
the SDSS data. The value ta = 11.2 kpc we find is the result 
of this compromise, but it is nevertheless unable to reconcile 
the model with the SDSS data. A similar compromise drives 
the fit in the estimates of log r\o ■ Because of the scaling with 
r~ 4 in the point mass potential, the term entering rn may 
contribute to the outer region circular velocity only if td is 
quite large. But, in such a case, R/rn will be much smaller 
than 1 in the inner regions thus making the contribution of 
this term to overcome by orders of magnitude the Newto- 
nian one in the regions where this latter is yet able to fit the 
data. As a consequence, r n must be small so that it does not 
contribute enough to boost v c with respect to the Newto- 
nian value hence motivating the need for a not too large ta ■ 
Indeed, we find an almost linear correlation between log tja 
and log rjD along the chain motivated by the fact that the 
larger is td, the larger is the contribution to v c (R) of the 
r~ 4 term of the potential and hence the smaller is the need 
for the r 2 term which translates in a larger ta- 

5 CONCLUSIONS 

Initially motivated by its attractive features from the point 
of view of quantum gravity, the HL proposal has soon be- 
come one of the most investigated theories of gravity. We 
have here complemented recent works on its cosmological 
consequences by addressing its impact on the gravitational 
potential. Contrary to the claim in TC09, we have demon- 
strated that static spherically symmetric solutions other 
than the Schwartzschild - de Sitter one exist. As a conse- 
quence, we have found a modified gravitational potential 
made out of the Newtonian 1/r one corrected by the addi- 
tion of three further terms scaling as r 2 , which corresponds 
to the effect of a cosmological constant, and a quickly de- 
creasing term, proportional to 1/r 4 . The importance of these 
two terms is parametrized in terms of two conveniently de- 
fined scaling radii, namely {ta, ro), related to the couplings 
entering the HL Lagrangian. In order to consider astrophysi- 
cally interesting situations, we have then developed a general 
formalism to compute the circular velocity curve provided 
the mass density and the mass function of the system are 
given. As an application, we have then evaluated the Milky 
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Way rotation curve using as only source of the gravitational 
field a spheroidal truncated power - law bulge and a dou- 
ble exponential disc. It turns out that the modified rotation 
curve is unable to fit the data thus demonstrating that the 
HL theory can not play the role of an alterative solution to 
the missing mass problem. 

It is worth noting that Mukohyama (2009) has shown 
that classical solutions to the infrared (IR) limit of the HL 
theory can mimic General Relativity plus cold dark matter 
so that one can argue that it should be possible for the HL 
theory to provide an effective dark matter halo in apparent 
contradiction with our finding. Actually, this is not the case. 
First, Mukohyama only refers to the HL IR limit so that the 
Lagrangian he has considered does not include any poten- 
tial term, while we here explicitly include this term through 
the coupling parameters gi. Second, the cold dark matter 
term comes out as a consequence of the global Hamiltonian 
constraint being less restrictive than the local one typically 
imposed in General Relativity. In our static spherically sym- 
metric metric, the global Hamiltonian constraint reduces to 
Eq.(5). In order to find our solution, we then impose that 
the integrand in Eq.(4) identically vanishes so that we are 
actually converting the global constraint in a local one thus 
going back to a situation similar to the General Relativ- 
ity case. As a consequence, we argue that the possibility to 
get a matter term as an integration constant is lost in our 
approach. 

The next step in our analysis of the HL theory on galac- 
tic scale should naturally be the inclusion of a dark halo. 
There are however some subtle issues making such a logi- 
cal step forward not so easy to address. First, needless to 
say, there are few hints on which the dark halo mass density 
profile should be. All the models used in literature are mo- 
tivated by the outcome of numerical simulations, such as, 
e.g., the popular NFW (Navarro et al. 1997) and Einasto 
(Einasto 1965; Cardone et al. 2005), or evidences from rota- 
tion curve fitting, such as the isothermal sphere (Binney & 
Tremaine 1987) and the Burkert (Burkert 1995; Burkert & 
Salucci 2000; Borriello & Salucci 2001) profile. All the pre- 
vious works implicitly assume the validity of the Newtonian 
potential so that they are no more valid in a modified frame- 
work as the one we are using here. As a consequence, we 
have therefore to explore a wide class of density profile able 
to mimic most of the models in literature to finally select the 
most empirically motivated one. Moreover, as Eq.(25) shows, 
we also need to know the dark matter mass function which is 
completely unknown. Should the dark matter be composed 
of pointlike particles all having the same mass, we can adopt 
a Dirac S leaving the mass of the particle as an unknown or 
setting it according to some particle physics model. As a 
final consequence, the fit to the rotation curve should de- 
termine the three HL scaling parameters (va, td) and Ndm 
halo parameters so that severe degeneracies among these 
Ndm quantities may take place. In order to reduce them, a 
better way should be to consider external galaxies rotation 
curves assuming that the MF is the same as the MW one. 
In such a way, we could take advantage of the data probing 
the full radial range and not only the outer disc as done here 
with the MW sample. Moreover, such a test can also probe 
the universality of the scaling radii (ta^d) thus providing 
a further mandatory test of the HL model. 

As a final remark, we want to stress that, should such 



an analysis be successful, one has still to address a different 
issue. Let us suppose that we have indeed well fitted the 
rotation curves of a large sample of spiral galaxies thus de- 
termining a halo model and the values of (va, ro). One could 
then use Eqs.(ll) and (18) to infer constraints on the HL 
coupling parameters (go, . . . , gs)- Although we have thus five 
constraints (the two radii plus the two conditions imposed in 
the derivation of the potential) so being unable to determine 
all the eight quantities (go, ■ ■ ■ , gs), one could nevertheless 
try to see whether the allowed region of the parameter space 
is consistent with what is inferred from cosmological analy- 
ses (Dutta & Saridakis 2009) or Solar Sytem tests (see e.g. 
Iorio & Ruggiero 2009 and references therein). It is worth 
noting that such an attempt should give a consistent picture 
of the HL framework on all physical scales thus allowing to 
draw a deeper insight into its viability. 
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